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METHOD AND APPARATUS FOR MEASURING BLOOD LIPOPROTEIN 
LEVELS BY NMR SPECTROSCOPY 

FIELD OF THE INVENTION 
The invention relates to the measurement of 
lipoprotein levels in blood plasma or blood serum and, more 
particularly, the levels of low-density lipoproteins (LDL) > 
5 high-density lipoproteins (HDL) , very low-density 
lipoproteins (VLDL) and subclasses thereof. These 
lipoproteins account for the vast majority of the 
cholesterol found in blood. 

BACKGROUND OF THE INVENTION 
10 The importance of accurately measuring cholesterol 

levels in blood is well known. The federal government, in 
combination with more than twenty health organizations, has 
launched an aggressive campaign, through the National 
Cholesterol Education Program, to convince physicians and 
15 the general population of the dangers of high cholesterol 
levels in the blood. All persons are urged to have their 
cholesterol levels checked, and specific treatments are 
recommended based on the precise measured cholesterol 
level. In addition, treatments are not based solely on the 
20 total cholesterol level, but instead, on the level of LDL 
chol sterol. LDL cholester 1 appears to be the major cause 
of clogged arteries, wh reas HDL cholesterol aids in 
removing cholest rol deposits. A separate, and m re 
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expensive test is requir d to determin the level of LDL 
cholesterol and it is usually not c nducted unless the 
measured total cholester 1 level is at the borderlin or 
high risk levels. 
5 The most common methods for measuring cholesterol 

levels are notoriously inaccurate and the standard practice 
is to repeat the measurement a number of times when high 
levels are detected on the first measurement. Inaccuracies 
of 5% or more have been found in nearly half of the 
10 measurements made by testing laboratories and 15% of the 
measurements were inaccurate by an amount greater than 10%. 
These inaccuracies are inherent in the current measurement 
methods which require considerable handling of the blood 
and certain presumptions about the ratios of its 
15 constituent parts. 

Direct quantization of lipoprotein cholesterol is 
usually achieved by enzymatic assay of the individual 
lipoproteins, which are separated by ultracentrifugation, 
electrophoresis, or selective precipitation. There is 
20 great variability among the available separation methods in 
terms of accuracy, convenience, and cost* Generally, the 
most accurate methods are those involving 
ultracentrifugation, but these are very time consuming and 
expensive and therefore not suitable for largescale 
25 population studies. The most widely used alternative is an 
indirect method introduced by W. T. Friedewald, R. I. Levy, 
and D. S. Fredrickson, Estimation of the Concentration of 
Low-Density Lipoprotein Cholesterol in Plasma, Without Use 
of the Preparative Ultracentrifuge, Clin. Chem. 18, 499-502 
30 (1972) . In this procedure, plasma triglyceride (TG) and 
total cholesterol (TC) are measured by enzymatic assay. To 
a separate aliquot of plasma is added one of several 
reagents which selectively precipitates VLDL and LDL. After 
removing the precipitate by centrifugation, the supernatant 
35 is assay d for cholesterol to pr vide a measure f HDL 
chol st rol (HDL-C) . An estimate of VLDL cholesterol 
(VLDL-C) is then made by dividing the plasma triglyceride 
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level by five. The LDL cholesterol (LDL-C) concentration 
is then calculated by differ nc : LDL-C = TC - (HDL-C + 
VLDL-C) . Although this meth d is relatively rapid and 
inexpensive, there are several steps where experimental 
5 error can be introduced, particularly in the precipitation 
step. In addition, the accuracy of the analysis depends on 
the assumption that VLDL-C can be reliably estimated as one 
fifth the concentration of plasma triglyceride. When 
fasting samples are used, this is generally true, but other 

10 formulas have also been suggested to give more accurate 
values as described by D. M. DeLong, E. R. DeLong, P. D. 
Wood, K. Lippel, and B. M. Rifkihd, A Comparison of Methods 
for the Estimation of Plasma Low- and Very Low-Density 
Lipoprotein Cholesterol, J. Am. Med. Assoc. 256, 2372-2377 

15 (1986). 

It has also been shown that the major lipoprotein 
constituents could be further subdivided into subclasses 
based on further refinement of particle densities. Krauss 
et al, J. Lipid Research 23, 97-104 (1982), Atger et al., 

20 Clinical Chemistry 37, 1149-1152 (1991). The distribution 
of these subclasses within a major lipoprotein group may in 
itself provide further insight into risk analysis of CHD. 
Stossel et al, JAMA 260, 1917-1921 (1988). However, 
previous methods of determining the distribution profile of 

25 subclasses have been time consuming and unable to determine 
a number of subclass concentrations simultaneously. 

SUMMARY OF THE INVENTION 
The present invention relates to a method and 
apparatus for measuring the lipoprotein constituents of 

30 blood using a nuclear magnetic resonance (NMR) technique. 
More specifically, the method and apparatus includes 
acquiring proton NMR data from a sample of blood plasma or 
serum, processing the acquired NMR data to produce a 
chemical shift spectrum, and deconvoluting the spectrum in 

35 terms of the spectra f subclasses of the major classes of 
lipoprotein, to give the concentrati n f each f the 
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lipoprotein constituents and the distribution of subclasses 
f the constituents. It has been discovered that the 
spectrum is accurately repr sented by a linear c mbination 
of the spectra of plasma constituents into which the blood 
5 can be fractionated. The major constituents are commonly 
classed according to density as VLDL, LDL, HDL, 
chylomicrons and protein. The NMR spectral properties of 
the subclasses of these classes have been found to be 
virtually invariant from person to person. Thus, any 
10 differences in the NMR spectra are due entirely to 
differences in the amplitudes of the subclass spectra, 
which, in turn, is due to the concentrations of the 
subclasses and therefore the constituents in the blood. 

A general object of the invention is to provide an 
15 accurate and reliable measurement of the lipoprotein 
constituents of blood. Since the observed spectrum of a 
whole plasma sample can be closely simulated by 
appropriately weighted sums of the NMR spectra of the 
subclasses of its constituent classes, it is possible to 
20 extract the concentrations of these constituents in a 
sample by calculating the weighting factors which give the 
best fit between the sample spectrum and the calculated 
spectrum. The handling and processing of the sample is 
relatively simple compared to prior methods and there is, 
25 therefore, less opportunity for error. Furthermore, by 
including chylomicrons as a constituent the fasting 
requirement of previous methods is no longer required. 

Another object of the invention is to provide a method 
for measuring the lipoprotein constituents of blood at an 
30 economical cost and on a mass basis. The preparation of 
the sample is a trivial task and the actual NMR measurement 
is carried out automatically by an NMR spectrometer in five 
minutes or less. The deconvolution calculations are also 
carried out automatically by a computer which prints out a 
35 report that indicates the concentrations of all of the 
lipoprotein subclasses. The sums of the subclass 
concentrations falling within a particular density range 



PCT/US92/06216 

5 

giv the concentration of the lip prot in class 
corr sponding to that d nsity range • Furth rm r , the 
distribution of subclasses f constituents is developed for 
each constituent simultaneously. 
5 Another object of the present invention is to provide 

a method which is independent of environmental variables 
for determining the concentrations of lipoproteins in 
blood. By aligning the subclass reference spectra and the 
sample spectra to a control peak, the line shape analysis 
10 using the deconvolution process is rendered independent of 
environmental variables such as temperature and sample 
composition. 

The foregoing and other objects and advantages of the 
invention will appear from the following description. In 

15 the description, reference is made to the accompanying 
drawings which form a part hereof, and in which there is 
shown by way of illustration a preferred embodiment of the 
invention. Such embodiment does not necessarily represent 
the full scope of the invention, however, and reference is 

20 made therefore to the claims herein for interpreting the 
scope of the invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 
Figure 1 is a graph showing the chemical shift spectra 
of a representative sample of lipoprotein constituent 
25 subclasses; 

Figure 2 is a graph showing the chemical shift spectra 
of a first plasma sample and its lipoprotein constituents; 

Figure 3 is a graph showing the chemical shift spectra 
of a different plasma sample and its corresponding 
30 lipoprotein constituents; 

Figure 4 is a graph showing the Agarose Gel Filtration 
Profile of the samples used in Figures 2 and 3; 

Figure 5 is a series of graphs of the subclass 
distributions of the samples used in Figures 2 and 3; 
35 Figur 6 is a flow chart of the program for performing 

the method of the pr sent invention; and 
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Figur 7 is a block diagram £ the apparatus employed 
t practice the present invent! n. 

GENERAL DESCRIPTION OF THE INVENTION 

The present invention relates to the calculation of 
5 concentrations of lipoprotein constituents of blood plasma* 
Specifically, the present invention determines lipoprotein 
concentrations using the deconvolution of proton NMR 
spectra of plasma in much the same manner as U.S. Patent 
No, 4 r 933,844, the specification of which is hereby 

10 incorporated herein by reference as if set out fully* 

X H NMR spectra of human "blood plasma contain two 
prominent peaks centered at approximately 1.2 and 0.8 ppm 
(relative to the chemical shift standard, TSP) . These 
peaks arise from methylene (CH 2 ) and methyl (CH 3 ) protons, 

15 respectively, of plasma lipids. Each of these peaks is 
very heterogeneous in nature, consisting of overlapping 
resonances from protons of the several chemically distinct 
classes of lipids present in plasma: triglycerides; 
cholesterol; cholesterol esters; and phospholipids. These 

20 lipids are packaged together into three major classes of 
lipoprotein particles, which differ in the proportions of 
lipids which they contain. These lipoprotein particles 
also differ in density from which their names are derived: 
very low density lipoprotein (VLDL) , low density 

25 lipoprotein (LDL) , and high density lipoprotein (HDL) . 
These major classes of lipoprotein constituents may be 
further subdivided into subclasses. A subclass of 
lipoprotein particles comprises particles which have common 
physical properties, such as density, which permit a 

30 subclass to be fractionated from other subclasses and that 
exhibits NMR properties which are distinct from other 
subclasses. The NMR properties of one subclass may be 
distinct in a number of ways such as chemical shift or 
lineshape variations which make the subclass 

35 distinguishable from other subclasses. Subclasses 
distinguished upon density may be considered as a subclass 
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of the class of lipoprotein which contains particles of th 
subclasses density. 

Only that fraction of the lipids in thes lipoprotein 
particles that are in a fluid, mobile state (as opposed to 
5 an ordered liquid-crystalline state) contribute to the 
plasma lipid NMR resonances. The heterogeneity of these 
plasma signals is reflected by their complex lineshapes, 
which vary from person to person owing to variations of the 
plasma concentrations of the different lipoprotein 

10 particles, each of which has its own characteristically 
different NMR spectral properties. 

The method of the present invention allows the 
concentrations of lipoprotein particles (VLDL, LDL, HDL, 
and chylomicrons) of a plasma sample to be extracted from 

15 its l H NMR spectrum by a computer analysis of the lineshapes 
of its methyl and methylene signals. Use of the methyl 
signal alone, however has been found to be preferable. The 
method exploits the finding that this region of the 
observed plasma spectrum is accurately represented by a 

20 simple linear combination of the spectra of subclasses of 
the five major lipoprotein classes into which plasma can be 
fractionated by differential flotation ultracentrifugation. 
The five classes are differentiated on the basis of their 
density (in kg/L) and include: VLDL (density < 1.006); LDL 

25 (density = 1.006 to 1.063); HDL (density = 1.063 to 1.21); 
"Protein" (density > 1.21) and chylomicrons (density < 
0.940). The "Protein" constituent is the mostly 
protein-containing bottom fraction left behind after 
flotation of the lipoproteins. The inclusion of the 

30 chylomicron constituent eliminates the need for the blood 
sample to be taken from a fasting donor. 

The NMR spectral properties of these classes have been 
found to be quite similar from person to person. This is 
illustrated in Table 1 which is the result of a study 

35 conducted at th University f Wisconsin-Milwaukee and the 
Medical College of Wise nsin. 
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TABLE 1 

500 MHz NMR Parameters £ the Separated 
Lipoprotein Constituents of Plasma 

Parameter Mean +/- SD 

5 VLDL (n » 117) 

CH 2 Chemical Shift (ppm) 1.233 +/" 0.002 

CH 3 Chemical Shift (ppm) 0.839 +/- 0.002 

CH 2 Linewidth (Hz) 20.8 +/" 1-9 

CH a Linewidth (Hz) 16.3 +/- 0.8 

10 CHa/CHs Intensity Ratio 3.76 +/- 0.29 

LDL (n = 66) 

CH 2 Chemical Shift (ppm) 1.219 +/- 0.005 

CH 3 Chemical Shift (ppm) 0.822 +/- 0.002 

CH 2 Linewidth (Hz) 34.0 +/" 2.9 

15 CH 3 Linewidth (Hz) 21-1.+/- 1.6 

CH2/CH3 Intensity Ratio 1.27 +/~ 0.13 

HDL (n * 70) 

CH 2 Chemical Shift (ppm) 1.186 +/- 0.004 

CH 3 Chemical Shift (ppm) 0.796 +/" 0.003 

20 CH 2 Linewidth (Hz) 34.4 +/- 2.9 

CH 3 Linewidth (Hz) 20.0 +/- 0.8 

CHa/CHa Intensity Ratio 1.58 +/- 0.13 

PROTEIN (n = 111) 

CH 2 /CH 3 Intensity Ratio 0.37 +/- 0.10 

25 The small person-to-person variations in the 

lineshapes of the lipoprotein classes are caused by the 
subclass heterogeneity known to exist within each of these 
lipoprotein classes. Figure 1 shows the lineshapes and 
chemical shifts (positions) for a number of subclasses of 

30 lipoproteins. As shown in Figure 1, the chemical shifts 
and lineshape differences between the subclasses are much 
smaller than those between the major lipoprotein classes, 
but are completely reproducible. Thus, differences among 
the NMR signals from the plasma of individuals are caused 

35 by differences in the amplitudes of the lipid resonances 
from the subclasses present in the plasma, which in turn 
are proportional to their concentrations in the plasma. 

This is illustrated in Figures 2 and 3 in which the 
NMR chemical shift spectra of two substantially diff rent 

40 bl od plasma samples are shown. The spectral peak pr duced 
by methyl (CH 3 ) protons 60 (shown as a solid line) is shown 
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for the blood samples in Figures 2 and 3. The spectral 
peak 61 (shown as a dotted line) in Figures 2 and 3 is 
produced by th arithmetic sum of the NMR signals produced 
by the lipoprotein subclasses of the major classes VLDL, 
5 LDL, HDL, proteins and chylomicrons, as illustratively 
shown in Figure 1. It can be seen that the lineshape of 
the whole plasma spectrum is dependent on the relative 
amounts of the lipoprotein subclasses whose amplitudes 
change dramatically with their relative concentrations in 
10 the plasma sample. It is the invariant lineshape of the 
NMR spectra of the subclasses of plasma lipoprotein 
constituents across the entire population and the fact that 
these lineshapes may be arithmetically added to produce the 
lineshape of the blood plasma sample, which is the basis 
15 for the present invention. 

Since the observed CH 3 lineshapes of whole plasma 
samples are closely simulated by the appropriately weighted 
sum of lipid signals of its constituent subclasses of 
lipoprotein classes, it is possible to extract the 
20 concentrations of these constituents present in any sample. 
This is accomplished by calculating the weighting factors 
which give the best fit between observed blood plasma NMR 
spectra and the calculated blood plasma spectra. The 
process of NMR lipoprotein analysis is thus comprised of 
25 the following steps: (1) acquisition of an NMR "reference" 
spectrum for each of the pure constituent lipoprotein 
subclasses of plasma, (2) acquisition of whole plasma NMR 
spectra using measurement conditions identical to those 
used to obtain the reference spectra, and (3) computer 
30 deconvolution of the plasma NMR spectra in terms of the 
constituent subclasses to give the concentration of each 
lipoprotein constituent expressed as a multiple of the 
concentration of the corresponding lipoprotein reference. 
The plasma lineshape analysis is accomplished by 
35 calculating weighting coefficients f r each f the 
reference NMR spectra which minimize the sum of squared 
deviations between the observed plasma NMR spectrum and 
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that which is calculated by summing the weighted reference 
spectra. 

The inclusion of the subclasses f the major 
lipoprotein classes decreases the error between the 
5 calculated lineshape and the NMR lineshape, thus increasing 
the accuracy of the measurement while allowing for 
simultaneous determination of the subclass profile of each 
class. Because the differences in subclass lineshapes and 
chemical shifts are small it is important to correctly 

10 align the reference spectrum of each subclass with the 
plasma spectrum. The alignment of these spectra is 
accomplished by the alignment of control peaks in the 
spectra which are known to respond in the same manner to 
environmental variables, such as temperature and sample 

15 composition, as do the lipoprotein spectra. One such 
suitable alignment peak is the peak produced by CaEDTA, 
although other EDTA peaks may be utilized. By alignment of 
the spectra, the small variations in the subclasses 
lineshapes and chemical shifts may be exploited to produce 

20 higher accuracy and subclass profiles. 

DETAILED DESCRIPTION OF THE DRAWINGS 
Figure 1 is a graph of the chemical shift spectra of 
representative subclasses of plasma that serve as reference 
spectra in the plasma lineshape analysis. As shown in 

25 Figure 1, the spectra labeled V2 and V6 are of chylomicrons 
and VLDL respectively; the spectra labeled L2 (1.006 < 
density < 1.035) and L5 (1.035 < density < 1.063) are of 
constituent subclasses of the LDL major class of plasma; 
and the spectra labeled H2 (1.063 < density < 1.125) and H5 

30 (1.125 < density < 1.210) are constituent subclasses of the 
HDL major class of plasma. The subclasses shown in Figure 
1 are representative and further refinement of densities or 
sizes within a constituent lipoprotein class results in 
additional subclasses which may be incorporated as 

35 reference spectra. As shown in Figure 1, the reference 
spectra f subclasses within a lipopr tein class exhibit 
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substantial similarity to oth r reference spectra within 
that lipoprotein class. Furthenn re, the spacing between 
spectra within a class is closer than between spectra of a 
different class. Because of the very close proximity 
5 between reference spectra for subclasses within a 
lipoprotein class, proper alignment of the reference 
spectra to the sample spectrum is essential to prevent 
misregistration of subclass spectra. 

Figures 2 and 3 are graphs of the chemical shift 
10 spectra of two different blood samples and the 
corresponding subclass constituents obtained using the 
present invention. The lineshapes shown in Figures 2 and 
3 are for the methyl peak of plasma. As shown in Figures 
2 and 3, the calculated lineshape 61 (dashed line) of the 
15 methyl peak and the experimental lineshape 60 (solid line) 
closely correspond to one another. The subclasses utilized 
to form the calculated lineshape are also shown in Figures 
2 and 3. The protein component of the lineshapes is also 
shown in Figures 2 and 3. Figures 2 and 3 illustrate the 
20 variation between individuals of the distribution of 
constituent subclasses within a major lipoprotein class. 
For example, in Figure 3, the reference spectra for the 
subclasses LI, V6, H2 and H3 were utilized in calculating 
the lineshape whereas in Figure 2 none of these subclasses 
25 were required. Another example is the large V4 component 
present in Figure 2 which is not present in Figure 3. 

The absence of a VLDL component in the sample 
reflected in Figure 3 is also reflected in Figure 4, which 
is an Agarose Gel Filtration Profile showing the 
30 lipoprotein constituent concentrations for the samples 
reflected in Figures 2 and 3. The gel filtration process, 
the results of which are shown in Figure 4, separates the 
lipoproteins of the plasma by particle size. Because the 
size of the lipoprotein particles is related to their 
35 density, Figure 4 also indicates the relative densities of 
the lipoprotein constituents f the plasma. As shown in 
Figure 4 , the tracing 30 corresponds to the sample utilized 
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in Figure 2 and the tracing 31 correso nds to the sample 
utilized in Figure 3. As can be seen in Figure 4, the 
profile 30 has a large p ak corresponding to the presence 
of VLDL in the sample. The profile 31 has little if any 
5 peak corresponding to the presence of VLDL. These 
differences are reflected in Figures 2 and 3 by the 
presence of the large V4 subclass in Figure 2 and the 
relatively small V6 subclass in Figure 3. Figure 4 further 
illustrates the effect of different subclass distributions 
10 within a major lipoprotein class. As shown in Figure 4, 
the LDL peak of tracing 30 occurs at higher elution volumes 
than does the LDL peak for tracing 31. The higher elution 
volumes of the LDL peak indicates a higher concentration of 
smaller LDL subclasses in the sample represented in Figure 
15 2, as indicated by the presence of L4 and L6, than was 
found in the sample represented in Figure 3, which 
comprised L4 and LI components. 

Figure 5 is a series of graphs illustrating the size 
distributions within the major lipoprotein classes for the 
20 samples reflected in Figures 2, 3 and 4. The distributions 
were obtained utilizing an average of duplicate analyses of 
blood samples from two individuals. As shown in Figure 5, 
variations in the distribution of the constituent 
subclasses of the major lipoprotein classes results in the 
25 variations of the major constituent peaks as seen in Figure 
4. 

Figure e is a flow chart of the program which may be 
executed on a suitable computer for performing the method 
of the present invention. The program implements a 

30 nonnegative linear least squares regression with simplex 
optimization. The implementation of the program, as shown 
in Figure 6, is similar to that described in U.S. Patent 
No. 4,933,844, column 9, line 11 through column 11 line 15 
is incorporated herein by reference as if set out fully. 

35 However, the program as d scribed in Figure 6, limits the 
coefficients t positive values while only calculating the 
real values of the coefficients, thereby reducing the 
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matrices t h their original size. The program described 
in Figure 6, further adds the simplex optimizati n to 
further refine the coefficient values. As sh wn in Figure 
6, the program is capable of accepting a plurality of 
5 reference spectra representing subclasses of lipoprotein 
classes and protein. 

In the preferred embodiment, the NMR measurements are 
conducted at '250 MHz using an unmodified commercial 
spectrometer, model WM250 manufactured by Bruker 
10 Instruments, Inc. A fixed- frequency 5 mm l H probe is 
installed and the temperature controller is set to 23 *C. 
( +/" 0.5 'C). Field homogeneity is optimized by shimming 
on a sample of 99.8% D20 until the spectral linewidth of 
the HDO NMR signal is less than 0.6 Hz. The 90* RF 
15 excitation pulse width is set to a value of 5.5 +/- 0.2 
microseconds for the D20 measurement. 

Referring particularly to Figure 7, the spectrometer 
indicated by dashed line 10 is controlled by a digital 
computer 11. The computer 11 is sold under the trade name 
20 "ASPECT 2000" and it has a 24-bit word length and storage 
for 80 K words. It is particularly well suited for 
performing fast Fourier transformations and includes for 
this purpose a hard-wired sine table and hardwired multiply 
and divide circuit. It also includes a data link 12 to an 
25 external personal computer 13, and a direct-memory-access 
channel 14 which connects to a hard disc unit 15. 

The digital computer 11 also includes a set of 
analog-to-digital converters, digital-to-analog converters 
and slow device I/O ports which connect through a pulse 
30 control and interface circuit 16 to the operating elements 
of the spectrometer. These elements include an RF 
transmitter 17 which produces an RF excitation pulse of the 
duration, frequency and magnitude directed by the digital 
computer 11, and an RF power amplifier 18 which amplifies 
35 the pulse and couples it to the RF transmit coil 19 that 
surrounds sample tube 20. Th NMR signal pr duced by the 
excited sample in the presenc of a 5.875 Tesla polarizing 
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magnetic field produced by superconducting magnet 21 is 
received by a coil 22 and applied t an RF receiver 23. 
The amplified and filtered NMR signal is dem dulated at 24 
and the resulting quadrature signals are applied to the 
5 interface circuit 16 where they are digitized and input 
through the digital computer 11 to a file in the disc 
storage 15. 

After the NMR data is acquired from the sample in the 
tube 20, it is processed by the computer 11 to produce 
10 another file which is stored in the disc storage 15. This 
second file is a digital representation of the chemical 
shift spectrum and it is subsequently read out to the 
personal computer 13 for storage in its disc storage 25. 
Under the direction of a program stored in its memory r the 
15 personal computer 13 processes the chemical shift spectrum 
in accordance with the teachings of the present invention 
to print a report which is output to a printer 26. 

It should be apparent to those skilled in the art that 
the functions performed by the personal computer 13 and its 
20 separate disc storage 25 may also be incorporated into the 
functions performed by the spectrometer's digital computer 
11. In such case, the printer 26 is connected* directly to 
the digital computer 11. Prior to their measurement, the 
0.5 ml reference samples are removed from the refrigerator 
25 and allowed to rise to a temperature of 23* C. for a period 
of from ten minutes to two hours. A sealed coaxial insert 
(Wilmad, Cat.#WGS-8BL) containing an external standard used 
for field-frequency lock and normalization of the plasma 
signal amplitudes is placed into each plasma NMR sample 
30 tube before the spectrum is run. The composition of this 
standard insert is 0.008M TSP (sodium 3-trimethyl 
[2,2,3,3- 2 H 4 ] propionate), 0.6 mM MnS0 4 , 99.8% D 2 0. The D 2 0 
provides the field-frequency lock signal and the integrated 
area of the TSP resonance is used to normalize the 
35 amplitudes of the plasma lipid resonances to correct for 
variations in spectrometer detection sensitivity. The 
solution is doped with Mn 2+ to paramagnetically broaden the 
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normally sharp TSP resonance to make Its integrated area 
insensitive to small differences in field horn geneity and 
t shorten its Tl relaxation time to a value comparable to 
those of the plasma lipid resonances (200 to 500 
5 milliseconds) . The reference sample containing the coaxial 
insert is placed at a defined depth in the sample tube and 
placed in the spectrometer. The sample is spun at a rate 
of 20 Hz. After locking on the D 2 0 signal from the coaxial 
insert, a brief shimming of the z and z* gradient controls 
10 is performed using the NMR signal of the plasma water. 

The reference spectra is then acquired using a 
standard one-pulse sequence preceded by a one second 
selective decoupler presaturation pulse of the strong H 2 0 
resonance. A spatially selective composite 90* observation 
15 pulse (9O,-9O y -9O.,-90. y ) is used to minimize water 
suppression artifacts as described by A. Bax, "A Spatially 
Selective Composite 90o Radiofrequency Pulse", in J. Magn. 
R&son. 65, 142-145 (1985), although a normal 90- pulse also 
gives satisfactory results. The following acquisition 
20 parameters are used: 240 transients (4 dummy scans), 4 K 
data size, quadrature detection, 2800 Hz spectral width 
(9.9 to - 1.2 ppm), 0.73 sec. acquisition time, 1.0 sec. 
decoupler presaturation pulse (0.2 watt) at the H20 
frequency, 22 microsecond composite 90o pulse, and constant 
25 receiver gain for all spectra. The time-domain spectra 
(FIDs) of the four lipoprotein reference samples are 
digitized and stored on computer disk. 

The reference sample FIDs are processed identically to 
give the frequency-domain spectra used for the plasma 
30 lineshape fitting analysis. The processing operations of 
Fourier transformation, phasing, and baseline correction 
are accomplished using the standard commercial software of 
the NMR spectrometer (Bruker "DISNMR" program) . The FIDs 
are Fourier transformed using 16K data points after 
35 application of a 1.0 Hz linebroadening exp nential 
multiplication function. All spectra are scaled 
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id ntically. The spectra are then phase corrected t give 
pure absorption mode signal. 

The system is now ready to measure plasma samples. 
The procedure is virtually the same as that described above 
5 for measurement of the s reference samples. The same NMR 
spectrometer is used and it is set up to operate in the 
identical fashion used to acquire the lipoprotein reference 
spectra. The time domain spectrum (FID) of the plasma 
sample is acquired in the identical fashion as the 
10 reference spectra and it is processed in the same manner to 
produce a digitized representation of the blood plasma 
sample spectrum in the disk of the personal computer. The 
whole plasma spectrum is then accurately referenced to the 
sharp NMR resonance peak produced by the calcium complex of 
15 EDTA which is present in the sample. The sample spectrum 
and the reference spectra are shifted as needed to align 
the CaEDTA peak at 2.519 ppm on the horizontal scale. 

The mathematics used in the lineshape fitting process 
(i.e. non-negative linear least squares fit with simplex 
20 optimization of an unknown function in terms of a weighted 
sum of known functions} is well known and is described in 
many textbooks of numerical analysis and in articles such 
as p. J. Leggett, Numerical Analysis of Multicomponent 
Spectra, Analytical Chemistry 49 , 276-281 C1977) • A 
25 program for performing this function on a PC-AT computer is 
described by the flow chart in Figure 6. 

EXAMPLE 1 

Blood is collected from healthy subjects after a 12 to 
14 hour fast into Vacutainer Tubes (Beet on Dickinson, 

30 Rutherford, NJ) containing EDTA (final EDTA concentration, 
1 g/L) . Plasma is separated within 2 hours by 
centrifugation (2000 x g, 20 minutes) and stored at 4 # C. 
Plasma and lipoprotein lipids are analyzed by automated 
procedures at a hospital clinical laborat ry. Total 

35 cholesterol and triglyceride concentrations are measured 
enzymatically with a Hitachi 717 analyzer (Boehringer 
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Mannheim Diagnostics, Indianapolis, IN) . HDL— C is measured 
with an Ektachem 700 analyzer (Eastman Kodak, R Chester, 
NY) in the supernate obtained after precipitati n fa 
plasma aliquot with dextran sulfate (M, 50, 000) -Mg 1 * . 
5 Fasting plasma samples are fractionated into their 

lipoprotein subclass components according to density by 
sequential flotation ultracentrifugation at 4*C as 
described by Schumaker and Puppione, Methods in Enzymology 
128, 1-68 (1988). The following components are isolated: 

10 VLDL (d<1.006 kg/L), large LDL (d=l. 006-1. 035 kg/L), small 
LDL (d=l. 035-1. 063 kg/L), HDI* (d=l. 063-1. 125, HDL, 
(d=l. 125-1.21 kg/L), and Protein (d>1.21 kg/L). 
Chylomicrons (d<0.940 kg/L) are isolated from plasma 
samples of subjects fed a fat-rich meal according to the 

15 procedure of Hatch and Lees, Adv. Lipid Res. 6, 1-68 
(1968). 

The above lipoprotein components from several subjects 
are combined to provide the standard samples used to 
generate the reference spectra employed in the computer 

20 lineshape analysis of the plasma spectra. To ensure a 
uniform ionic composition, which is essential for correct 
alignment of the reference spectra, each lipoprotein 
component solution is dialyzed for 24 hours at 4°C against 
three changes of dialysate. The dialysate contains 120 

25 mmol KC1, 5 mmol of EDTA, 1 mmol of CaCl 2 , and 1 g of NaN 3 , 
pH 7.4) Each component is then concentrated at 4'C to 
about fivefold its normal plasma concentration using a 
Centricon-10 microconcentrator (Amicon, Inc. ) and then 
stored at 4'C prior to NMR analysis. 

30 All spectra of the isolated lipoprotein components and 

the real and artificial plasma samples are acquired under 
identical conditions at 250 MHz with a Bruker WM-250 
spectrometer (Bruker Instruments, Billerica, MA) . Samples 
(0.5 mL) in 5-mm (o.d.) NMR tubes are stored at 4°C for as 

35 long as six days before being analyzed. A systematic study 
of the effect of sample storage conditions on plasma methyl 
and methylene lineshape indicates that spectral changes of 
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samples kept at 4°C are negligible for the first six days, 
but occasionally are apparent after longer storag (notably 
for plasma with high concentrations of triglyceride) . 
Before placing samples in the spectrometer, they are 
5 allowed to equilibrate for 15-30 minutes at room 
temperature. A sealed coaxial insert, containing an 
external intensity standard (sodium 3- 
trimethylsilyl [2,2,3, 3-fy] propionate) , was used as 
previously described in Otvos et al., Clin. Chem. 37, 369- 

10 376 (1991) . Each NMR sample is placed at a reproducible, 
defined depth in the proton probe and allowed to 
equilibrate for 5 minutes at the chosen sample temperature 
(15 - 45 °C). The probe is detuned by several megahertz to 
prevent radiation damping, which increases the 90° pulse 

15 length from 6 to 16jis. Spectra are run locked with the 
sample spinning (20Hz) and the magnetic field homogeneity 
is optimized for each sample by shimming on the water 
signal. A spatially selective composite 90° observation 
pulse is used to minimize water suppression artifacts, 

20 although a normal 90° pulse also gives satisfactory 
results. The spectral width is set to 2800 Hz, the data 
size is 4K r the acquisition time is 0.73 sec, the 
composite pulse length is 64 fis, and the number of 
transients is 120 with 4 dummy scans and constant receiver 

25 gain. The time-domain data are zero-filled to 16K 
multiplied by a 1-Hz exponential linebroadening function, 
and Fourier-transformed with identical scaling. After 
phasing and chemical shift referencing to the sharp CaEDTA 
resonance at 2.519 ppm, a linear baseline was applied as a 

30 correction to flatten the baseline between 1.8 and -0.2 
ppm. 

The Fourier-transformed plasma spectra and those of 
the pure lipoprotein reference samples are stored on the 
magnetic disk of an IBM-compatible PC-AT computer after 
35 transfer from the Aspect 2000A computer of the Bruker 
spectrom t r by use of the FASTRAN program (University of 
Wisconsin-Madison) . The linear least-squares analysis of 
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the lineshape f the plasma methyl lipid res nance is 
performed by using a program written in BASIC. This 
program first places the real and imaginary data p ints 
from the methyl region of the plasma and lipoprotein 
5 reference spectra into separate arrays in computer memory. 

Several additional "approximated 1 ' reference spectra 
are also added into computer memory to account for known 
lipoprotein subclass size/density heterogeneity (and hence 
spectral heterogeneity) beyond that already accounted for 
10 by the subclasses isolated to provide the lipoprotein 
reference spectra (chylomicrons, VLDL, large and small LDL, 
HDI*, HDL 3/ and protein). The "approximated" reference 
spectra are created by digitally shifting the methyl 
regions of the lipoprotein component spectra to the left or 
15 right (downfield or upfield) by an appropriate amount. 
Thus, the chylomicron spectrum (V2) is shifted two data 
points (0.68 Hz) to the left and right to create artificial 
spectra VI and V3, respectively, to represent a larger and 
smaller population of chylomicron particles. Similarly, 
20 the VLDL spectrum (V6) is shifted to the left by two data 
points (V4) and one data point (V5) and to the right by one 
data point (V7) to represent a wider range of VLDL particle 
sizes. The large LDL spectrum (L2) is shifted left (LI) 
and right (L3) by 1 data point as is small LDL (L5) to give 
25 spectra representative of "larger" (L4) and "smaller" (L6) 
small LDL. The HDI* spectrum (H2) was shifted left (HI) and 
right (H3) by two data points and the HDLj spectrum (H5) one 
data point to the left (H4) and right (H6) to give a total 
of six HDL subclass spectra representative of the range of 
30 HDL particle sizes expected. Thus, including the spectrum 
of the Protein component (d>1.21 kg/L) , the plasma 
lineshapes are fit using a total of 20 reference spectra: 
V1-V7 representative of chylomicron and VLDL constituents, 
L1-L6 representative of LDL subclasses, and H1-H6 
35 representative of the HDL subclasses (in the numbering 
system used, th larger numbers designate subclasses of 
smaller particle diameter r greater density) . 
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The lineshape deconvolution is achieved with a 
nonnegative linear least-squares program described by 
Lawson et al., Solving Least Squares Problems , (Prentice 
Hall ,1974) that introduces the physical constraint that the 
5 derived concentrations must be positive. The latter 
constraint is necessary when fitting plasma samples that 
may not contain one or more of the components included in 
the fit, because experimental errors in the data (noise) 
often cause the calculation to give negative concentrations 

10 for these components. For example, including a chylomicron 
component in the analysis of fasting plasma samples will 
frequently give negative chylomicron concentrations (and 
hence incorrect concentrations for the other lipoprotein 
constituents) when the unconstrained least-squares method 

15 is used. 

Mathematically, the methyl lineshape analysis is 
described by the following equation: 

where the superscripts R and I denote the real and 
imaginary parts of the spectra; P, is the experimental 

20 plasma spectrum, consisting in this case of 132 discrete 
data points; are the reference spectra of the n 

lipoprotein components; V H is the spectrum of the "protein 19 
component; and c ir c*, and c £ are the unknown relative 
concentrations whose values are determined by minimizing 

25 the root mean square deviation between the experimental 
plasma spectrum and the calculated spectrum. 

The relative lipoprotein concentrations, c i# derived by 
this method have no absolute meaning since they only relate 
the concentrations of the lipoprotein components of the 

30 plasma sample to those of the reference components of 
arbitrary concentration. How ver, if the concentrations of 
the lipopr tein solutions used to provide the standard 
referenc spectra have been accurately determined by 
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chemical analysis (ie., by chol sterol analysis or, in the 
case of VLDL, by triglyceride analysis) the relative 
concentrations, c i§ derived from the lineshap analysis can 
be readily multiplied by these chemical concentrations to 
5 give lipoprotein concentrations expressed in the usual 
terms of lipoprotein cholesterol or triglyceride 
concentrations. If total chylomycron, VLDL, LDL, and HDL 
concentrations are desired, they are obtained simply by 
adding the derived concentrations of the individual 
10 subclasses (V1-V3, V4-V7, L1-L6, and H1-H6 respectively). 
Note that, if chylomicrons are treated as a subclass of 
VLDL, then VLDL concentration is obtained from VI -V7. If 
profiles of the lipoprotein subclass distributions are 
desired, they are provided directly by the relative 
15 concentrations of the subclass components used in the 
lineshape fitting algorithm. Alternatively, a single "size 
distribution parameter" that gives the weighted average 
particle size within a given class of lipoproteins may 
readily be calculated. 
20 The methyl resonance lineshape of chylomicrons is very 

similar to that of VLDL, but the signal is shifted slightly 
downfield. To determine the effect of the presence of high 
concentrations of chylomicrons on the quantification of 
VLDL, LDL, and HDL by lineshape fitting, we analyzed by NMR 
25 and chemical methods both fasting and nonfasting plasma 
samples from individuals fed a fat-rich meal. The presence 
of chylomicrons in the postprandial sample is easily 
discerned by the altered position of the plasma peak 
maximum. By including chylomicrons as a fifth component in 
30 the lineshape analysis, where only four other components 
corresponding to the major lipoprotein classes were used, 
we obtained an excellent fit of the experimental plasma 
spectrum, but also found the derived concentrations of 
VLDL, LDL, and HDL to be nearly identical to those in the 
35 fasting state (see Table 3 below). These and similar 
results obtained for other postprandial samples indicate 
that lipoproteins can be reliably analyzed in nonfasting 



22 



PCT/US92/06216 



TG 


TO 


HDL-C 


LDL-C 


2.21 


1.88 


0.41 


1.03 


4.47 


1.91 


0.40 


0.66 


6.20 


1.87 


0.37 


026 



plasma by NMR. In contrast , the widely used method f 
Fried wald et al., Clin. Chem. 18, 499-502 (1972) has an 
absolute r quirement for fasting samples because the 
accuracy of LDL-C values is severely compromised by the 
5 presence of chylomicrons. 

Table 2 

Influence of Chylomicrons on the NMR 
Lipoprotein Assay 

Plasma Lipid Concentration. o/L 

10 Sample 
Fasting 

2 h postprandial 
4 h postprandial 

NMR-derived lipoprotein concn. mmol/L proton 
15 Sample Chylomicrons VLDL LDL HDL 

Fasting — 19.0 las 10.6 

2 h postprandial 22.4 19.2 13.3 9.9 

4 h postprandial 32.4 17.9 1Z9 10.8 



The information derived from the above procedure, 

20 which is very rapid (minutes) and requires almost no sample 
manipulation, is equivalent to that provided by acquiring 
separate spectra of the components prepared by 
ultracentrifugation (days) and comparing the integrals of 
their lipid NMR signals to those of reference lipoprotein 

25 samples. It is important to note that what is being 
measured by this procedure (NMR signal amplitude 
originating from the "mobile" lipid molecules in each class 
of lipoprotein) is related to, but fundamentally different 
from, lipoprotein lipid and protein concentrations derived 

30 by the various chemical and immunochemical assays in 
current clinical use. There is thus no reason to expect a 
perfect correlation to exist between these NMR-derived 
lipoprotein levels and those derived from standard serum 
cholesterol and triglyceride analyses. Despite well 

35 document d limitations in the accuracy and precision of the 
latter measurements, they are in widespread clinical use 
because of their proven value in assessing coronary heart 
dis ase risk and other lipid-related dis ase states. It is 
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possible that lipoprotein levels derived from the NMR 
lineshape deconvolution proc ss may have even greater 
diagnostic utility, but this will n t be kn wn until 
extensive clinical correlation studies have been performed. 
5 It should be apparent to those skilled in the art that 

many variations are possible from the above-described 
preferred embodiment of the invention. For example , the 
polarizing field strength may be increased to further 
spread the NMR spectrum and to thereby improve the 

10 resolution of the deconvolution process. Also, the 
measurements may be conducted at other temperatures. 
Regardless of the magnetic field strength or the 
measurement temperature which is chosen, it is important 
that the chosen values remain constant throughout the 

15 process of producing the reference spectra and the sample 
spectra. 



PCT/US92/06216 

24 

That which is claimed is: 

1. A method of measuring the lipoprotein 
constituents of blood, comprising: 

storing the NMR spectra of a plurality of lipoprotein 
5 classes as reference spectra for said classes, said storing 
step further comprising storing the NMR spectra of a 
plurality of subclasses for at least one of said 
lipoprotein classes as reference spectra for the subclasses 
of said lipoprotein class; 
10 acquiring an NMR signal produced by a plasma or serum 

sample in an NMR spectrometer; 

producing an NMR spectrum of the sample by 
transforming the acquired signal; 

producing a calculated lineshape by adding together 
15 the stored reference spectra in amounts determined by 
respective reference spectrum coefficients; and 

adjusting the reference coefficients to fit the 
calculated lineshape to the NMR spectrum of the sample. 

2. A method of measuring the lipoprotein 
20 constituents of blood according to Claim 1 further 

comprising the step of calculating the concentration of at 
least one major lipoprotein class or subclass thereof as a 
function of the value of the reference coefficients. 

3. A method of measuring the lipoprotein 
25 constituents of blood according to Claim 1 in which the 

lipoprotein subclasses are selected from the group 
consisting of the subclasses of VLDL, LDL, HDL and 
chylomicrons. 

4. A method of measuring the lipoprotein 
30 constituents of blood according to Claim 1 in which the NMR 

spectrum includes the peak produced by methyl protons. 
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5. A method of measuring the lip pr tein 
constituents of blood according t Claim 1 in which the 
calculated lineshape is fit to the NMR spectrum of the 
sample by minimizing the root mean square error. 

5 6. A method of measuring the lipoprotein 

constituents of blood according to Claim 1 in which the 
calculated lineshape is fit to the NMR spectrum of the 
sample through nonnegative linear least squares 
deconvolut ion . 

10 7. A method of measuring the lipoprotein 

constituents of blood according to Claim 1 further 
comprising the step of aligning the position of the sample 
spectrum to the position of the reference spectra.' 

8. A method of measuring the lipoprotein components 
15 of blood according to Claim 7 wherein said alignment step 

comprises aligning a control peak in the reference spectra 
to the corresponding control peak of the sample spectra, 
wherein said control peak is a non-1 ipoprotein constituent 
peak which responds to environmental variables in the same 
20 manner as lipoproteins. 

9. A method of measuring the lipoprotein components 
of blood according to Claim 8 wherein said control peak is 
an EDTA peak. 

10. A method of measuring the lipoprotein 
25 constituents of blood, comprising: 

storing the NMR spectrum of lipoprotein constituents 
selected from the group consisting of VLDL, LDL, HDL, 
chylomicrons and protein as a reference spectrum for that 
constituent, said storing step further comprising storing 
30 the NMR spectra of a plurality of subclasses of at least 
one of the lipoprotein constituents of plasma as reference 
spectra for the subclasses of said lipoprotein constituent; 
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acquiring an NMR signal produc d by a plasma or serum 
sample in an NMR spectrometer; 

producing an NMR spectrum of the sample by 
transforming the acquired signal; 
5 aligning the position of the sample spectrum to the 

position of the reference spectra; 

producing a calculated lineshape by adding together 
the stored reference spectra in amounts determined by 
respective reference spectra coefficients; and 
10 adjusting the reference coefficients to fit the 

calculated lineshape to the NMR spectrum of the sample. 

11. A method of measuring the lipoprotein 
constituents of blood according to Claim 10 further 
comprising the step of calculating the concentration of at 

15 least one major lipoprotein class or subclass thereof as a 
function of the value of the reference coefficients. 

12. A method of measuring the lipoprotein 
constituents of blood according to Claim 10 in which the 
NMR spectrum includes the peak produced by methyl protons. 

20 13. A method of measuring the lipoprotein 

constituents of blood according to Claim 10 in which the 
calculated lineshape is fit to the NMR spectrum of the 
sample by minimizing the root mean square error. 

14. A method of measuring the lipoprotein components 
25 of blood according to Claim 10 in which the calculated 

lineshape is fit to the NMR spectrum of the sample through 
nonnegative linear least squares deconvolution. 

15. A method of measuring the lipoprotein components 
of blood according to Claim 10 wherein said alignment step 

30 comprises aligning a control peak in the reference spectra 
to the corresponding control p ak of the sample spectrum, 
wherein said control p ak is a non-lipoprotein constituent 



27 



PCT/US92/06216 



peak which responds to environmental variables in the same 
manner as lipoproteins. 

16. A method of measuring the lipoprotein components 
of blood according to Claim 10 wherein said control peak is 

5 an EDTA peak. 

17. A method of measuring the lipoprotein 
constituents of blood, comprising: 

storing the NMR spectra of a plurality of lipoprotein 
constituents as reference spectra for said constituents; 
10 said storing step further comprising storing the NMR 
spectra of chylomicrons as reference spectra for 
chylomicrons; 

acquiring an NMR signal produced by a plasma or serum 
sample in an NMR spectrometer; 
15 producing an NMR spectrum of the sample by 

transforming the acquired signal; 

producing a calculated lineshape by adding together 
the stored reference spectra in amounts determined by 
respective reference spectrum coefficients; and 
20 adjusting the reference coefficients to fit the 

calculated lineshape to the NMR spectrum of the sample. 

18. A method of measuring the lipoprotein 
constituents of blood according to Claim 17 further 
comprising the step of calculating the concentration of at 

25 least one major lipoprotein constituent as a function of 
the value of the reference coefficients. 

19. A method of measuring the lipoprotein 
constituents of blood according to Claim 17 in which the 
NMR spectrum includes the peak produced by methyl protons. 

30 20. A method of measuring the lipoprotein components 

of blood according to Claim 17 in which the calculated 
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lineshape is fit to the NMR spectra of the sample through 
nonnegative linear least squares deconvoluti n. 

21. A method of measuring the lipoprotein 
constituents of blood according to Claim 17 further 

5 comprising the step of aligning the position of the sample 
spectrum to the position of the reference spectra; 

22. A method of measuring the lipoprotein components 
of blood according to Claim 21 wherein said alignment step 
comprises aligning a control peak of the reference spectra 

10 to the corresponding control peak of the sample spectrum, 
wherein said control peak is a non-1 ipoprotein constituent 
peak which responds to environmental variables in the same 
manner as lipoproteins. 

23. A method of measuring the lipoprotein components 
15 of blood according to Claim 22 wherein said control peak is 

an EDTA peak. 

24. A method of measuring the lipoprotein 
constituents of blood, comprising: 

storing the NMR spectra of a plurality of lipoprotein 
20 classes as reference spectra for said classes; 

acquiring an NMR signal produced by a plasma or serum 
sample in an NMR spectrometer; 

producing an NMR spectrum of the sample by 
transforming the acquired signal; 
25 aligning the position of the sample spectrum to the 

position of the reference spectra by aligning a control 
peak in the reference spectra to the corresponding control 
peak of the sample spectra, wherein said control peak is a 
non-lipoprotein constituent peak which responds to 
30 environmental variables in the same manner as lipoproteins; 

producing a calculated lineshape by adding together 
the stored reference spectra in amounts determined by 
respective referenc spectrum coefficients; and 
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adjusting the r ference coefficients t fit the 
calculated lineshape to th NMR sp ctrum of the sample. 

25. An apparatus for measuring the lipoprotein 
constituents of blood/ comprising: 
5 means for storing the NMR spectra of a plurality of 

lipoprotein classes as reference spectra for said classes, 
said storing step further comprising storing the NMR 
spectra of a plurality of subclasses for at least one of 
said lipoprotein classes as reference spectra for the 
10 subclasses of said lipoprotein class; 

means for acquiring an NMR signal produced by a plasma 
or serum sample in an NMR spectrometer; 

means for producing an NMR spectrum of the sample by 
transforming the acquired signal; 
15 means for producing a calculated lineshape by adding 

together the stored reference spectra in amounts determined 
by respective reference spectrum coefficients; and 

means for adjusting the reference coefficients to fit 
the calculated lineshape to the NMR spectrum of the sample. 

20 26. An apparatus for measuring the lipoprotein 

constituents of blood according to Claim 25 further 
comprising means for calculating the concentration of at 
least one major lipoprotein class or subclass thereof as a 
function of the value of the reference coefficients. 

25 27. An apparatus for measuring the lipoprotein 

constituents of blood according to Claim 25 further 
comprising means for aligning the position of the sample 
spectrum to the position of the reference spectra. 
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(START) 



n=# OF. LIPOPROTEIN COMPONENTS 
m'=n+1, m"=n+2 
m=n+3 
q=# OF DATA POINTS 
rms-10 6 
rmss=10 6 
CVG=10~6 (OR SMALLER) 

* 

READ REAL (R) AND IMAGINARY 
(I) PLASMA SPECTRAL DATA INTO 
P* AND vRj, i=1 TO q 

t 

READ LIPOPROTEIN COMPONENT 
SPECTRA INTO VR , 
j=1 TO n, i=1 TO Jl q 

* 

READ REAL AND IMAGINARY 
PROTEIN COMPONENT SPECTRAL- 
DATA INTO VR/; AND V^/Zj, 
i=1 TO q 

SET FLAGS N(j) TO INCLUDE OR 
EXCLUDE COMPONENTS IN FIT: 
N(j)=1 THEN INCLUDE Vj 
N(j)-0 THEN EXCLUDE Vj 
FIX cj =0 

* 

STORE ORIGINAL N(j) IN NS(j) 
1 = # OF ns(j) NOT EQUAL TO 0 
(USED WITH SIMPLEX) 

FIG. 6A. 
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CALL SOLVE SUBROUTINE 

♦ • ■ ■ 



CALL RMSD SUBROUTINE 





CD- 


NO 


chsq>10 i:) 
? 










Tyes 






SET- 


N(k)=0 
.ck=0. 


-FOR ALL c k <0 














CALL SOLVE SUBROUTINE 














CALL RMSD SUBROUTINE | 


YES I 


chsq>10 ID 
« — . ? 1 









[NO 

FORM m+1 VERTEX VECTORS 
WITH ORIGINAL c" VECTOR: 
Sk = c, k=1 TO m+1 

AC|< = 0, k=1 TO m 
IF NS(k)=1 THEN 
Ac k =.01, k=1 TO m 

FIG. 6B. 
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1 



VARY ith PARAMETER TO CREATE 
m+1 VERTICES IN m SPACE FOR 
SIMPLEX » Sjj=sjj+Aq, i=1 TO m 



I 



k=0 



rsd(k)=chsq 

I " 



k=k+1 



CALL RMSD 
SUBROUTINE 



m 

p?«5>jk Vji 


NO 


j=i 




i=1 TO q 






DETERMINE VERTICES WITH 
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